The Nature of Explosive Percolation Phase Transition 
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In this Letter, we show that the explosive percolation is a novel continuous phase transition. 
The order-parameter-distribution histogram at the percolation threshold is studied in Erdos-Renyi 
networks, scale-free networks, and square lattice. In finite system, two well-defined Gaussian-like 
peaks coexist, and the valley between the two peaks is suppressed with the system size increasing. 
This finite-size effect always appears in typical first-order phase transition. However, both of the 
two peaks shift to zero point in a power law manner, which indicates the explosive percolation 
is continuous in the thermodynamic limit. The nature of explosive percolation in all the three 
structures is belong to this novel continuous phase transition. Various scaling exponents concerning 
the order-parameter-distribution are obtained. 
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I. INTRODUCTION 

Percolation [l[ , the simplest model presenting continu- 
ous phase transition, is one of the fundamental problems 
in statistical physics, since it provides deeper understand- 
ing of many other issues through Fortuin-Kasteleyn rep- 
resentation . The percolation model itself has applica- 
tions to a wide variety of different systems, ranging from 
sol-gel transition and polymerization [!, to conduc- 
tivity of composite materials and flow through porous 
media la Jab to spreading process and robustness in net- 
works [7h12|. Hitherto the critical properties in most of 
these systems are well described by the universality of 
percolation model in corresponding dimensionality. 

Strikingly, Achlioptas, D'Souza, and Spencer |l3[ re- 
ported that the percolation transition for the Erdos- 
Renyi (ER) model |14j may become discontinuous, 
through a modified growth procedure known as prod- 
uct rule (PR). They found at the percolation threshold 
an abrupt jump in the size of the largest component, 
which was named as explosive percolation (EP) compared 
with the traditional continuous percolation transition. In 
light of this, subsequent studies were devoted to uncover- 
ing the underlying mechanism of EP [l5j , proposing new 
models for EP [lq421| . and studying EP transition with 
different topologies and dimensionalities |22H26| . Re- 
cently, two empirical studies focused on the EP in human 
protein network [23] and social network [28| . 

While further investigations confirmed the abrupt 
transition in EP, it was also shown that the critical dis- 
tribution of cluster sizes follows a power law 26| . which 
manifests the features characteristic of the second-order 
phase transition. The contradictions make the nature of 
explosive percolation transition a controversial issue in 
statistical physics, which needs to be clarified. Recent 
research proves that the explosive percolation is a weak 
continuous phase transition in mean-field structure [29r - 
l32l ]. It is reminiscent of the weak first-order phase tran- 
sition [33| in five-state Potts model [34| , where, since the 
correlation length is very large at transition point, the 



accessible system size in numerical simulation is always 
in the critical region, and thus the picture of cluster dis- 
tribution is characterized by fractal shapes rather than 
smooth droplets. The weak first-order phase transition 
is hard to establish due to its proximity and resemblance 
to a critical point, and weak continuous phase transition 
is hard to confirmed due to the smallness of the critical 
exponent for the order parameter. 

However, in this Letter, we show that the explosive 
percolation is not a traditional weak continuous phase 
transition by examining the distribution histogram of 
the order parameter G defined as the fraction of ver- 
tices in the largest cluster. Three structures are consid- 
ered, such as ER network, scale-free (SF) network, and 
two-dimensional (2D) lattice. Three key features are ob- 
served in all these structures. Firstly, we find that in 
finite system two well-defined Gaussian-like peaks coex- 
ist in the order-parameter-distribution histogram at the 
percolation threshold, which represent the nonpercola- 
tive phase and percolative phase, respectively. Secondly, 
the probability of realizing a configuration in the inter- 
mediate phase between the two peaks is suppressed as a 
power law with the system size increasing. Finally, both 
of the two peaks shift to zero point in power law man- 
ner. These observations indicate the explosive percola- 
tion is a continuous phase transition with first-order-like 
finite-size effect. Various scaling exponents concerning 
order-parameter-distribution are obtained. 



II. METHOD 

For concreteness, numerical simulations were per- 
formed according to the original PR process [l3| : In each 
turn, two unoccupied edges are randomly chosen; the one 
which minimizes the product of the masses of the clusters 
it joins is retained. For square lattice, we imposed peri- 
odic boundary conditions in both directions to reduce the 
boundary effect. For SF network, we adopted the model 
by Chung and Lu (CL) j36j to build the network. Specif- 
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FIG. 1. Plots of A(G,p) as a function of order parameter G 
in the vicinity of percolation threshold p c for PR model (a) 
(b) (c) and traditional RG model (d) (e) (f). The exponent 
A = 2.8 is used for SF network. The system size is 4096 for all 
the simulations in this figure. The full information of p can 
be found in the movies of order parameter distribution in the 
vicinity of transition point [35l |. 



ically, every vertex in the system is assigned a weight 
beforehand according to the desired degree distribution, 
and at every time step, two edges are independently se- 
lected with probability proportional to the product of the 
weights of the vertices at the end of each edge. Then the 
PR is used to decide which is the next occupied edge. 

The controlling parameter p denotes the number of 
added edges divided by the system size N. We measured 
the order-parameter-distribution histogram H(G,p) for 
each p through extensive Monte Carlo (MC) simulations. 
According to the standard probability theory, the number 
of realized configurations with order parameter G is 

H{G,p) = exp[-A(G,p)} ~ Z- 1 (p)Q(G,p), (1) 

where Z(jp) is the normalization factor and Q(G,p) is 
the order-parameter probability density function, i.e., the 
probability that, after pN edges are added with PR pro- 
cess, the fraction of vertices in the largest cluster is G. 
When the number of realizations increases to infinity, 



2.5 
2.0 
1.5 

A 

1.0 
0.5 
0.0 









2.0 r 


- i 


(a) ER network 




1.8 - 




//\\ 




1.6 - 




#1 




AA 








1.4 - 




1 li/nL u 

1 /// NT 1 1 




1.2 - 



n 2.0 








h 1-6 
\ M 






V 1.2 






\ 1C -9 -8 -7 




n -InW 


(d) ER network 




0.0000 0.0004 - 0.0008 

N 




> lattice 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.0000 0.0003 ,0.0006 

G AT 1 




0.0000 0.0004 . 0.0008 

AT 1 



FIG. 2. The dependence of A(G,p) at percolation thresh- 
old on the order parameter G for different system sizes (left 
panels). The values of the minima have been normalized 
to 0, and the curves have been translated along horizon- 
tal axis for better comparisons. The system sizes N = 
1024, 2048, 4096, 8192, 16384 for (a) ER network and (c) SF 
network, and TV = 1600, 2500, 4900, 8100, 16900 for (b) 2D lat- 
tice. Right panels display the depth of the minima in A(G,p c ) 
as a function of system size. The red lines in the insets are 
the linear fits. The exponent A = 2.8 is used for SF network. 



H(G,p) is identical to Q(G,p) multiplied by a constant. 
Intuitively, we have A(G,p) — —\nH(G,p), and thus at 
a given p the location of the global minimum in A(G,p) 
denotes the most probable size of the giant component. 



III. RESULT AND DISCUSSION 

Our computer implementation makes use of the effec- 
tive Newman- Ziff algorithm [37} for tracking the largest 
cluster in the system. We carried out 10 6 MC sweeps per 
vertex to achieve high statistical accuracy for H(G,p). 
In Fig. [TJ we show the behavior of A(G,p) near the per- 
colation threshold p c for both PR model and traditional 
random growth (RG) model. It is well known that the 
percolation transition with RG is continuous, which is 
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FIG. 4. The locations of the two minima G~ (left) and G 
(right) as functions of system size N for scale-free networks 
and ER network. 



FIG. 3. The inset shows in SF network the size-dependent 
depth of minima in A(G,p) at percolation threshold for dif- 
ferent A. The red solid squares represent the same relation 
in ER network. The fitting exponent 9 as a function of A is 
displayed in the main panel. 



reproduced in the simulations (see Fig. [T](d) ( e ) if))- As 
p passes through the percolation threshold, there is only 
one global minimum in A(G,p), which implies the or- 
der parameter grows continuously from one phase to the 
other. For PR model, the situation is completely differ- 
ent. As p goes below the critical value, a local minimum 
appears in the region of large order parameter, and its 
value gradually approaches that of the global one. Right 
at the percolation threshold p c , the two minima have 
equal depth, indicating that the nonpercolative and per- 
colative configurations are realized with equal probabil- 
ity. When p is beyond p c , the second minimum becomes 
global and percolative phase dominates. The physical 
picture of the whole process [35| is reminiscent of the 
Landau theory of first-order phase transition. This is the 
first finite-size property of explosive percolation. 

In the following, we perform detailed finite-size scaling 
analysis of A(G,p c ). An important quantity is the depth 
of the minima (corresponding to the peaks in H(G,P)), 
AA, relative to the local maximum (corresponding to 
the valley in H(G, P)) in-between corresponding to the 
intermediate phase. For typical first-order phase transi- 
tions |38|, with the system size increasing, A A also 
monotonically increases as the minima gradually develop, 
and eventually goes to infinity in the thermodynamic 
limit. We show here that this property also exists in ex- 
plosive percolation. Since p denotes the number of added 
edges which is not continuous, it is more convenient to 
calculate the depth with A A = A max - (A 1 min + A 2 mm ) /2, 
where A max is the value of the local maximum and j4^ in 
and 4^ m are those of the two minima. For the determi- 
nation of A] nin and 4^ in , it should be sufficient to use 
quadratic fit in the vicinity of the minima corresponding 
to the gaussian-like peaks in H(G,p c ). However, A(G,p c ) 
shows sizable asymmetry near the minima, thus we use 
cubic fit to obtain more accurate results. To determine 



the value of A max , we found it adequate to use the same 
fit. Figure [2] shows the simulation results for the size- 
dependent behavior of A(G,p c ). Indeed, the depth of the 
minima monotonically increases with system size, and 
tends towards infinity in the thermodynamic limit. Fur- 
thermore, a clear relation between AA and the logarithm 
of N is observed, 



AA 



nnN. 



(2) 



In other words, the relative probability of finding a con- 
figuration in the intermediate phase is suppressed, as the 
system size increases, in power-law manner with expo- 
nent 9. It should be noticed that 9 E r = 0.310(1) for ER 
network and O2D = 0.097(1) for square lattice are very 
different. For SF network, this scaling relation holds for 
A > 2.0 (see the inset of Fig. [3]). As A increases, the value 
of exponent 9 gradually approaches that for ER network, 
and at A — > 00 the CL model is identical to ER model. 
In fact, 9(X) is already saturated for A > 3.0 (see Fig. [3]), 
where with PR the SF network generated by CL model is 
hardly distinguishable from ER network [24], [26| . The ex- 
ponent 9 for different structures is listed in Table Q] This 
is the second finite-size property of explosive percolation. 

The locations of the two minima in order-parameter di- 
mension, G + and G~ , directly measures the jump of the 
order parameter at percolation threshold. In all the three 
structures, both of G + and G~ shift to zero in a power 
law manner, G ~ rj^ 1 . As N increases, the two minima 
gradually get close to each other, and in the thermody- 
namic limit they merge at the transition point of the 
order parameter, where it presents a continuous phase 
transition. Actually, in the research of explosive percola- 
tion on scale-free network where another field, the degree 
distribution exponent A, comes into play, it was claimed 
that there exists a tricritical point (TP) at A c 6(2.3, 2.4), 
above which the EP transition is first order 24 1. At the 
same time, however, careful finite-size scaling analysis 
implied that for A < 3.0 the EP transition is continu- 
ous [25[. Our results indicate that explosive percolation 
in scale-free network is continuous in the whole range 
of A. The scaling behavior of G ± in scale- free networks 
are displayed in Fig. 01 As A increases, the value of ex- 
ponents rj^ gradually approaches those for ER network. 
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TABLE I. Summary of scaling exponents obtained from our 
numerical analysis. 
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SF A = 2.6 


0.351(2) 


0.106(1) 


0.302(1) 


SF A = 2.8 


0.343(2) 


0.091(1) 


0.306(1) 


SF A = 3.0 


0.331(2) 


0.079(1) 


0.309(1) 


SF A = 3.2 


0.293(1) 


0.047(1) 


0.311(1) 


SF A = 3.5 


0.291(2) 


0.041(1) 


0.311(1) 


SF A = 4.0 


0.290(2) 


0.038(1) 


0.311(1) 


ER network 


0.282(1) 


0.035(1) 


0.310(1) 


2D lattice 


0.085(3) 


0.013(1) 


0.097(1) 



The exponents ry for different structures are listed in 
Table H 

IV. CONCLUSION 

In this Letter, we show that the explosive percolation is 
a novel continuous phase transition in ER network, scale- 



free network, and 2D lattice. By examining the order- 
parameter- distribution histogram at percolation thresh- 
old, it is found that two well-defined Gaussian-like peaks 
coexist, which represent the nonpercolative phase and 
percolative phase, respectively. Moreover, the probabil- 
ity of realizing a configuration in the intermediate phase 
between the two peaks is suppressed as a power law with 
the system size increasing. On the other hand, two peaks 
gradually get close to each other with the system size in- 
creasing, and in the thermodynamic limit they merge at 
the transition point of the order parameter. These obser- 
vations suggest the explosive percolation is a continuous 
phase transition with first-order-like finite-size effect. 
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